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(57) Abstract: A method of recognising a radar target comprises producing a sequence of Doppler spectra of radar returns form a 
scene and producing therefinom a sequence of Doppler feature vectors for a target in the scene. Hidden Markov modelling (HMM) 
is then used to identify the sequence of Doppler feature vectors as indicating a member of a particular class of targets. HMM is 
used to identify die sequence of Doppler feature vectors by assigning to each feature vector an occurrence probability by selecting a 
probability distribution or state from a set Uiereof associated with a class of targets, multiplying the occurrence probabilities together 
to obtain an overall probability, repeating for other probability distributions in the set to determine a combination of probability 
distributions giving highest ovaall probability for that class of target, then repeating for at least one other class of targets and selecting 
the target class as being that which yields the highest overall occurrence probability. 
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METHOD AND APPARATUS FOR RECOGNISING A RADAR TARGET 

This invention relates to a method of and an apparatus for recognising a radar target 

5 l\4oving target identification (MTl) radars are known. Such radars can provide all-weather, day 
and night surveillance capability and have a v\dde-area search capability as a result of rapid 
scanning of the radar beam. An MTl radar may be used to provide location information 
relating to a moving target and may also be used to identify a target as belonging to a 
particular dass of targets. For example an MTl radar may be used to classify a target as 

10 being a person, a wheeled vehicle or a traclced vehicle. Presently, the task of target 
recognition using an MTl radar is carried out by a human operator. In one known MTl radar, 
echo signals from targets are converted to audio signals which are output to an audio speaker. 
To identify a particular target, the radar's antenna is directed to a target for several seconds 
during which time an operator listens to the audio speaker output. An operator requires 

15 .; training to develop the ability to classify targets from their corresponding audio frequencies. 
However, if a scene under observation contains many different types of targets, a human 
operator, however.trained. cannot classify all types in the scene and cannot provide up-tondate 
information on a rapidly changing scene. 

20 In the case of imaging radars, to perfomi target recognition on radar signals algorithms have 
been developed which rely on the availability of 2-dimensional high-rBSolution imagery of 
targets. They exploit the differences between spatial structures of radar images in ortfer to 
classify targets. They are unfortunately unsuftable for use with MTl radars which are mainly 
low resoluticin sensors: MTl radars show target images as single bright point objects with no 

25 spatial structure in either range or azimuth. 

There are certain circumstances In which an MTl radar can achieve better spatial resolution. If 
the radar operates at higher bandwidth then the radar signature will have a higher slant range 
resolution and the target will be resolved in range, e.g. a bandwidth of 500 MHz will give a 
30 slant range resolution of 0.3 m. However, It is not possible to get higher azimuth resolution. 
This is because imaging radars use a synthetic aperture radar (SAR) processing, which 
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provides an azimuth resolution that is many times smaller than the real beam azimuth 
aperture. 

SAR processing assumes a static scene and moving objects will not be focused in azimuth 
using it Thus at best an MTI radar can only obtain 1 -dimensional high range resolution 
profiles of moving targets. Limited success has been achieved in developing recognition 
algorithms that exploit the 1 -dimensional high range resolution profile of targets for 
classification. A technique based on 1 -dimensional range profile template matching is 
described in the Proceedings of the International Radar Symposium IRS-98. There are, 
however, two main drawbacks with a range profile classifier. Firstly, high bandwidth radars 
are more expensive and it is more difficult to extract high range resolution profiles because 
moving targets can migrate through a succession of range cells and causing considerable 
difficulties. Secondly, target range profiles are very sensitive to changes In target orientation 
to the radar and the target physical build characteristics. 

Since a target can have any number of articulations, a range of build variations and an almost 
Infinite number of external fitting configurations, the number of potential Independent 
realisations of range profiles of a particular target Is large. This presents a major problem for 
classifiers using template matching, because they carry out recognition by comparing the 
range profile of an unknown target with a set of templates and choosing the class that yields 
the best match: this relies on having a template set that contain every single Independent 
realisation of target signature for a class of targets to be recognised. 

Another option open to an MTI radar is to use the Doppler characteristics of moving taitfets in 
determining the target classification. Doppler Is the phenomenon by which the radar return 
from an object is shifted In frequency due to the object's radial motion relative to a rad^r 
system. A Doppler return ft-om a target can be observed by looking at a received echo in the 
frequency domain. Raw radar data firom a target is recorded as a series of temporal samples. 
Using a n-point Fast Fourier Transfomi (FFT), temporal samples fl^om a target are transformed 
into a spectmm comprising n frequency samples or bins. The n frequency samples form a 
Doppler profile where the maximum unambiguous fl^quency is given by the Inverse of the 
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radar pulse repetition interval and tiie Doppler resolution per frequency bin is 1/nth of the 
unambiguous Doppler. 

A target classifier might be based on template matching similar to a range profile classifier but 
using templates of independent realisations of tiie Doppler profiles of moving targets. 
Altiiough Doppler profiles are comparatively less sensitive to target articulations, ttiey still 
exhibit fluctuations with changing target orientation and imaging geometiy. Furfliermore, 
Doppler profiles vary as a function of target speed and a profile's shape is modulated as a 
result of a target's vibrational and rotational motion. A target classifier based on template 
matching would require a prohibrtively large number of stored Doppler profiles. 

The present invention provides a metiiod of recognising a radar target comprising receiving 
radar returns fl-om a scene and processing the returns to produce a Doppler spectrum, 
characterised in that the metiiod also includes' processing tiie Doppler specbxim to produce a 
Doppler feature vector and using hidden Markov modelling (HMM) to identify the Doppler 
feature vector as indicating a member of a particular class of targets. 

The invention provides ttie advantage that moving targets may be classified wittiout recourse 
to a large dataliase of radar signature data. 

In a preferred embodiment, the method includes arranging for targets to be encompassed 
witiiin a single radar range cell. It may involve processing radar retums to obtain a sequence 
of Doppler spectra for each target and producing therefrom a sequence of Doppler feature 
vectors, and using HMM to Identify ttie seque.nce of Doppler feature vectors as indicating a 
member of a particular class of targets. This enables classification to use linl<ed Doppler 
profiles as part of a structural observation sequence: it exploits th^ fact that radar data for a 
target provides a series of Doppler profiles each oflset slightiy in time. Each Doppler profile is 
difl^erent, but over a sequence of observations tiie shape of the profile varies accordingly tp 
some detemiinistic process. This embodiment exploits useful infbnmation in ttie variation of the 
Doppler profile witii time. HMM algorithms have evolved to make use of sequences of speech 
signals in tiiis way. 
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The method may include using HMM to identffy the sequence of Doppler feature vectors by 
assigning to each feature vector an occun-ence probability by selecting a probability 
distribution or state from a set thereof associated with a class of targets, multiplying .the 
occurrence probabilities together to obtain an overall probability, repeating for other probability 
distributions in the set to detemiine a combination of probability distributions giving highest 
overall probability for that dass of targets, then repeating for at least one other class of targets 
and selecting tiie target dass as being tiiat which )de)ds the highest overall occun^ence 
probability. Probability distributions for successive feature vectors, may be selected on tiie* 
basis of some transitions between distributions allocated to successive feature vectors being 
allowed for the class of targets and others being forbidden. 

A preliminary HMM training procedure may be implemented in which parameters for the states 
or probability distributions and ti^nsition probabilities are produced by deriving Doppler feature 
vectors for training data obtained from known dasses of targets and calculating the mean and 
variance of vectors con^esponding to like targets or target dasses. Ttie HMM training* 
procedure may indude a plurality of cycles through state sequences. 

In another aspect, the present invention provides an apparatus for recognising a radar target 
comprising a radar receiver for receiving radar returns from a scene and processing means for 
processing ttie returns to produce a Doppler spectium characterised in that tiie processing 
means also indudes means for processing the Doppler spedrum to produce a Doppler feature 
vector and hidden Markov modelling (HMM) means for identifying the Doppler.feature vector 
as indicating a member of a particular dass of targets. The apparatus may be arranged to 
encompass targets within a single range cell. 

The processing means may be arranged to process radar returns to obtain a sequence of 
Doppler spedra for each target and to produce therefrom a sequence of Doppler feature 
vectors, and the HMM means may be arranged to identify the sequence of Doppler feature 
vedors as indicating a member of a particular dass of targets. 

The HMM means may be arranged to identify .tiie sequence of Doppler feature vectors by 
assigning to each feature vedor an occun-ence probability by selecting a probability 
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distribution or state from a set thereof assodated with a class of targets." multiplying the 
occurrence probabilities together to obtain an overall probability, repeating for other probability 
distributions in the set to detennine a combination of probability distributions giving highest 
overall probability for that class of targets, then repeating for at least one other class of targets 
and selecting the target dass as being that which yields the highest overall occun-ence 
probability. 

The HIWM means may be arranged to select probability distributions for successive feature 
vectors on the basis of some transitions between distributions allocated to successive feature 
vectors being allowed for the dass of targets and others being forbidden. It may be an-anged 
to undergo a preliminary training procedure in which parameters for the stat^ or probability 
distributions and transition probabilities are produced by deriving Doppler feature vectors for 
training data obtained from known dasses of targets and calculating the mean and variance of 
vectors corresponding to like targets or target dasses. The HMM training procedure may 
indude a plurality of cydes through state sequences. 

In order that the Invention might be more fully understood, embodiments thereof vAW now be 
described, by way of example only, with reference to the accompanying drawings. In which: 

Figure 1 is a schematic block diagram of a radar target recognition system of ttie 
invention; 

Figure 2 is a flow diagram of processing routines used in the Figure 1 system; 

Figure 3 illustates use of an LDA routine to calculate an LDA ti-ansfonnation function in the 
Figure 1 system; 

Figure 4 illusti^ates modilying target files by LDA conversion in ttie Figure 1 system; 
Figure 5 illusti^tes HMM parameter determination; 
Figures 6 to 8 lllustiBte HMM topologies; and 

Figure 9 is a schematic block diagram of anoUier embodiment of aradar target recognition 
system of the invention. 
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Referring to Figure 1 there is shown a radar target recx)gnition system of the invention 
indicated generally by 10. The system 10 comprises an antenna 12, a radar device 14, an 
analogue-to-digital (A to D) converter 16, a data recorder 18. a computer 20 and a. display 
device 22. The computer 20 is under operator control via a mouse 24 and a keyboard 28. Data 
passes between parts of the system 10 along data lines 30, of which those between the 
computer and A-to-D converter 16. the computer 20 and data recorder 18, and the A-to-D 
converter 16 and data recorder 18 include respective switches 33, 34 and 35. Control lines 36 
provide for computer-generated control signals to pass to the radar device 14, the A-to-D 
converter 16 and the display device 22. A communication line 37 is provided for passing 
messages between the computer 20 and data recorder 1 8. 

The antenna 12 is a fixed, high gain dish antenna. The radar device 14 is a "Midas" radar 
which is manufactured by Thales-Wells and described in the paper "Wide bandwidth mobile 
radar for ISAR/SAR radar imaging" by J D Mckenzie and E J Brown-Kenyon in the 
Proceedings of the Institute of Electrical Engineers, volume 219 pages 1 to 6, 1994. 

The Midas radar is a static ground based coherent pulsed radar. It uses the same antenna for 
transmitting and receiving and is therefore a raonostatic system. It operates at a centre 
frequency of 15.75GHzand at this frequency the antenna has a circular beam width of 19 
degrees.* The radar has a peak transmit power of 200W that gives it a maximum operating 
range of 2000m between itself and the target. It is capable of operating at a bandwidth ranging 
from 10MHz to 500 MHz. For this embodiment the operating, bandwidth is set to 15.6MHz 
which results in a range resolution of 9.6m. The radar is linearly polarised with both transmit 
and receive polarisations set to horizontal. It operates by transmitting a series of pulses at an 
interval of 2.441X1 0"* seconds. 

After each pulse has been transmitted, the radar is quiescent for a short duration and then 
records 32 values conresponding to radar echo from a series of 32 concun-ent range cells. The' 
range cells are arranged along a radial line extending outwardly from the antenna centre. The 
first range cell is the nearest in range to the antenna and the last the furthest in range. The 
radar uses the time delay between the end of transmission of a pulse and the start of the 
recording of the first range cell to determine the slant range offset to the start of the first range 
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cell. Each range cell contains the radar backscatter from a region bounded by a non- 
overlapping footprint on tennln being monitored by the radar. Any object or part of an object 
and all parts of the ten^ln that are exposed to the radar and fall within the boundary of the 
range cell radar footprint contribute to the radar echo received by the radar. The radar foo^rint 
corresponding to each range cell can be approximated to a rectangular shape with a width 
along slant range equal to the radar range resolution, which Is 9.6m and a cross-range width 
given by the product RsjOaz: heree„ is the antenna azimuth beam width that is 1.9 degree in 
this embodiment and Rs, Is the slant range offeet from the radar to the ith range cell. Rs, is a 
function of the range offset to the first range cell and the gate number and is given by:- 

Rs,= To + Ar(i-1). 1^1^32 (1) 

where Ar is the slant range resolution and ro is the range offset to the first range cell. For this 
embodiment ro is set to 260m. The range offset for the 32™* range cell is therefore 557.2m. The 
radar footprint for each range cell has a constant range extent but increases linearly with slant 
range In the azimuth direction. For this embodiment the azimuth extent increases from 8.3m 
for the first range cell to 18.5m for the 32™* range cell. TTie maximum dimensibn of targets to 
be classified by this embodiment is 8m the radar is capable of placing the radar echo from the 
entire target into the value of a single range cell. It is not in fact essential to do this but it 
makes the recognition process simpler - the target can be spread over a number of rahge 
cells if desired. 

The antenna can pan 360 degrees in azimuth and from -10 degrees to +20 degrees in 
elevation. The elevation angle is set to 0 degrees and the azimuth angle is adjusted manually 
so that It Is pointing directly at the target. The radar will capture the backscatter from the target 
when It comes within the slant range interval of 557.2m to 250m from the radar. Each of the 
32 values that the radar records for each pulse that is transmitted is a complex valve with a 
real and imaginary part. The analogue signal is passed through a A-to-D converter where the 
signal is digitised. Ail subsequent processing is performed on the digital data. The radar has a 
pulse repetition interval is 2.441x10"^ seconds, and it transmits 4096 pulses per second and 
receives data for 32 range cells for each pulse. 
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The A-to-D converter 16 fs of a standard type, for example a Tektronics model TKADIOC and 
is capable of digitising both real and imaginary parts of a complex input signal. The data 
recorder 18 is a standard high-speed magnetic tape recorder, for example an Ampex 107 
DCRsI recorder which records data at a rate of 107 IVIbits s'\ The computer 20 is a standard 
parallel computer, for example a Hewlett Packard model V2500 having twenty CPU nodes of 
type PA8500, a dock rate of 440 MHz and a memory capacity of 0.5 terabytes. The system 
10 has a graphical user interface (GUI) which is displayed on the display device 22 and with 
which an operator may Interact with the system 10 using the mouse 24 and the keyboard 28, 
Results generated by the system 10 are also displayed on the display device 22, together vwth 
standard information generated by radar apparatus, such as a target's range and velocity. 

The system 10 recogniseis and classifies moving targets as belonging to one of three target 
groups, namely personnel, wheeled vehicles and tracked vehicles. The vehicles must be of 
spedfic shapes and In standard configurations, for example vehide doors must be dosed. A 
target is reliably recognised by the system 10 if it is moving in a straight line either head-on or 
tail-on with respect to the antenna 12. A target is said to move head-on to the antenna 12 if it 
is moving directly towards it, and is said to move tail-on to the antenna 12 if the target is 
moving directly away from the antenna 12. If the target is a person. rBcognition will reliably 
occur if the person is moving at walldng or jogging pace. If the target is a wheeled vehide or a 
tracked vehide, recognition will reliably occur if the vehide is travelling a speed of 20 km h'\ 

Before the system 10 can be used to recognise and dassify unknown targets moving in the 
manner described above into the three target groups, it must be trained to do so. The system 
1 0 is trained as follows. 

The system 10 is first used to gather training data conresponding to radar retums flrom a set of 
co-operative example targets which comprises target classes into which the system 10 is 
required to dassify unknown targets in practice. In the present embodiment, the target groups 
are personnel, wheeled vehides and tracked vehicles. Switches 33 and 34 are opened and 
switch 35 is closed. Members of a set of co-operative example targets comprising personnel, 
wheeled vehides and tracked vehicles are brought serially into the main radiation lobe of the 
antenna 12. The co-operative targets move head-on and tail-on with respect to the antenna 
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12. The personnel move at walking and jogging paces, and.the wheeled and tracked vehicles 
move at a speed of 20 km h \ Twenty-five co-operative example targets are used to train the 
system 10. Table 1 shows the number of co-operative example targets belonging to each 
target group and the numbers that are anranged to move head-on and taii-on to the antenna 
12. 



Target Group 


Head-on 


Tail-on 


Personnel 


6 


2 


Tracked Vehicles 


5 


4 


Wheeled Vehicles 


4 


4 



Table 1 



Training data gathered by the system 10 for a particular target group includes data from just 
one example target within that group. For example, for the target group "wheeled vehicles", 
training data may be gathered from a lony: However, in use, the system 10 is able to 
recognise and classify any wheeled vehicle for example, a car, a Land Rover or a lony. 

The co-operative tanget is set into motion and the radar antenna is steered to point directly at 
the target The radar is switched on and starts transmitting pulses at 4 kHz pulse repetition 
frequency. A radar operator monitors the received radar echo on display device 22: when the 
target comes into the 250m to 557.2m range interval, the operator switches on the radar 
record mode. The operator uses the keyboard 28 to enable the record mode. Once the 
record mode has been enabled, for each pulse that is transmitted the radar device 1'4 captures 
data for 32 range cells which are digitised by the A-to-D converter 16 and stored serially on 
the data recorder 18. 

One pulse of the radar generates 32 complex values that are stored as a continuous data 
block on the data recorder 18 which is sandwiched between a header block and a footer block. 
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The data for the header block and footer blod< is generated by radar device 14 and contain 
information regarding the radar operating mode set up and recording mode set up which are 
useful for post processing purposes. Data for a total of 16384 pulses is recorded after which 
the recording mode is disabled. The data recorder device 14 sends tape addresses for the 
start and end of the recording sequence to the computer 20 via the communicafion line 37 
which is displayed on the display device 22 using the data link 30. The operator takes a 
manual record of the start and end tape addresses and also notes the target type and whether 
it was moving head-on or taikon. .This Information Is required in order to form a target file. 

The process is repeated for a total of twenty-five co-operative targets that are listed in Table 1 
and that manoeuvre in the manner explained earlier: i.e. vehicles moving at 20 km/hr and 
personnel at walking or jogging pace. Once all the twenty-five co-operative targets have been 
imaged the data capture phase is complete. 

On completion of data capture, switch 34 is closed and switch 35 is opened. Using keyboard 
28 the operator enters a tape address for the first co-operative target. The computer 20 sends . 
an instruction to the data recorder 18 via the communication link 37 to read the data from the 
tape at the entered address and process it through a series of routines which run on the 
computer 20. . 

Referring now to Figure 2 there is shown a flow diagram Illustrating the series of routines 51 to 
54 executed by the computer 20 to process data from the data recorder 18. The series 
comprises a data decoding routine 51, a target extraction routine 52, a Doppler normalisation 
routine 53 and a temporal to Doppler frames conversion routine 54. At the entered address in 
the data recorder 18 the computer 20 reads dgta for 16384 pulses. The data format used to 
record data on the data recorder 18 consisted of a header block followed by a data block of 32 
complex values corresponding to the 32 range cells recorded per pulse followed by a footer 
block. The decode routine 51 strips of the header and footer infomiation and places the 32 
complex values into a data matrix. The decode operation is performed for all the 16384 
pulses that were recorded on the data recorder 18 during data capture. The output after the 
decode routine 51 is a 2-dimensional matrix that has 16384 rows with each row containing 32 
columns. Each row represents radar data received from individual pulses and each column 
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represents data from individual range cells. Therefore, for example element (100, 12) will be 
the radar data derived from the 1 00th pulse at the 1 2th range cell. 

The 2-dlmentional data matrix from decode routine 51 is processed by the target extraction 
routine 52, The "extract" or target extraction routine 52 detects the range cell that contains the 
radar retum from the co-operative example target This is achieved using standard detection 
algorithms based on fixed thresholding or constant false aianm raite adaptive thrBsiiolding 
which are familiar to those skilled in the art of radar engineering. Once the range cell occupied 
by the co-operative example target has been identified, the 2-dimensfonal data matrix is 
reduced in size by dropping alj columns with the exception of the target range cell column. 
This results in reducing the 2-dimenslonal data matrix to a l-dimenslonal anay that has 16384 
rows and Just one column. The 1-dimenslonal anay Is data from a single range cell occupied 
by the co-operative example target over 16384 pulses. The 16384 element an^y is an array 
of temporal samples of the co-operative example target where each element Is separated 
temporally by 2.441x1 0"* seconds. The total array is equivalent to a temporal duration of 
2.441x10-^ X 16384 = 4 seconds. 

The 1-dimensional array from routine 52 Is then passed through to the Doppler normalisation 
routine 53. The input to the routine 53 is a 16384 element temporal anray. A16384 point FFT 
is perfomned on the 1-dimensional data that produces a 16384 element Doppler spectrum. 
The 1-dimensipnal array now represents frequency response from the co-operative example 
target Each element in the anay Is a frequency or Doppler bin. The first element 
con-esponds to a frequen<y of 0 Hz and the last element 4 kHz. Intennediate bins are 
separated in frequency by 4096/16384 Hz. Since the co-operative example target has a non- 
zero radial velocity with respect to the radar antenna 12 it will manifest as a peal< in the 
Doppler spectmm somewhere between frequency bin numbers 0 and 16383. Using a peak 
detector the Doppler bin containing the peak amplitude is Identified. The whole 16384 array is 
then shifted to the I6ft by the peak amplitude Doppler bin number so that the peak in the 
Doppler spectnjm con-esponds to Doppler bin zero. This has the effect of normalising the data 
so that the target Doppler appears to be 0 Hz. An inverse FFT is perfomned on the resultant 
16384 array to retum to an array that represents temporal (time domain) samples of the co- 
operative example target. 
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The 16384 element temporal array from routine 53 is then processed by the data Conversion 
routine 54, transforms the 1-dimensional temporal anay into a 2-dimensional Doppler frame 
matrix. In routine 54, each of the 16384 elements is initially a complex number with real and 
imaginary components: the imaginary component is discarded leaving only the real 
component The first 128 elements in the 16384 element array are then, selected and 
transfomried Into a 128 element Doppler vector using a 128-polnt FFT with Hanriming weighting 
window. The 128 elements in the Doppler vector are complex values that are converted into 
log-amplitude values. The Doppler vector constitutes one Doppler frame of the radar signature 
of the co-operative example target The Doppler frame is placed into the first row of a new 2- 
dimenslonal matrix. The 126 element of the. Doppler frame form the elements of the 128 
columns of the 2-dimensional matrix. Then another 128 elements are selected from the 
original 16384 element anray of temporal samples using an offset of 20 elements from where 
the previous selection of 128 elements was made. These 128 temporal samples are 
processed In a similar manner to generate a Doppler Frame vector. This second frame is 
placed along row 2 of the new 2-dimensional matrix. This process is repeated, each time 
selecting another 128 element vector from the array of temporal samples using the 
appropriate off set, perfomning the FFT and storing the Doppler frame in the 2-dimensionaI 
matrix until there are no more samples to select from the original 16384 element temporal 
sample anray. Through the choice of the size of the FFT window and the overiap in the 
elements used In calculating the Doppler frame vectors the process produces 800 
independent Doppler frames. Thus the final size of the 2-dimensional matrix Is 800 rows by 
128 columns. Each row In this 2-dimensjonal matrix constitutes data for the co-operative 
example target for an Independent Doppler frame where the elements of the columns are the 
actual Doppler spectral values. 

The 2-dimensional matrix output by* convert routine 54 is stored as a target file 25 within the 
computer 20. Using the keyboard 28. the radar operator gives the target file 25 a file name: 
The naming convention is letter/underscore/letter/underscore/numben the first letter in the 
name represents target class, P for personnel, T for fraclced vehicles and W for wheeled 
vehicles; the second letter represents target direction of motion, H for head-on and T for tail- 
on; the number uniquely identifies each target file from a particular class. The numbers are 
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not repeated for the target files belonging to a particular class: they may however be reused 
when naming target files that belong to separate classes. Thus a target file could be named * 
P_ H_ 2, which denotes that the target file is from an example of a co-operative target of the 
type personnel that is moving head-on to antenna 12 and it is file number 2 in this particular 
class. Since each target file 25 contains a 2 dimensional matrix in which each row is an 
independent Doppler frame, the entire file therefore represents a sequence of Doppier frames 
from the co-operative example target. 

Target files for the other twenty-four co-operative example targets are generated in a lilce 
manner. A total of twenty-five target files are ttierefore stored in the computer 20, each 
representing a sequence of Doppler frames that constitutes the target signature of an example 
co-operative target 

Each target file 25 consists of a series of 800 Doppler frames (tiiey are the rows of the 2- 
dimensionai matrix) and each Doppler frame is a 128 element vector. However, the elements 
in the Doppler frame are not independent and the size of the Qoppler frame vector can be 
reduced without loss of inforniation. This embodiment of the invention uses an algorithm 
called Linear Discriminant Analysis (LDA) to perfomi tfie data reduction. IDA are a widely 
used aigorittim and a good description of tiiis can be found in a paper by M. J. Hunt and S.iVI. 
Richardson in 'An Investigation of PLP and IMELDA acoustic representations and their 
potential combination', in proceedings of tiie IEEE ICASSP 1991, volume 2, pp 881 to 884. 
LDA is meant to improve inforniation content by exaggerating differences between natural 
groupings in the data set (i.e. between classes) while inhibiting differences witfiin the groups. 
The LDA process calculates a linear transformation that converts an n^lement vector to a m- 
element vector where m<n. This can be expressed symbolically as 

Y = gG9 (2) 

Where X is the n-element vector, Y is the m-element vector and g() is ttie linear transfomnation 
function. This embodiment of the invention uses the 'Rsher Criteria' LDA which states that the 
maximum number of elements of tiie vector Y_is one less the total number of categories. The 
LDA process categories are shown in Table 2: 
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Table 2 



Person 


Head-on 


Person 


Tail-on 


Wheeled Vehicle 


Head-on 


Wheeled Vehicle 


Taik)ri 


Tracked Vehicle 


Head-on 


Tracked Vehicle 


Ta[I-t)n 



There are a total of six categories. These are termed flne! label since for recognition only 
three dass labels are defined namely, 'personnel', tracked vehicle and 'wheeled vehicle'. 
Since the invention uses 'Fisher Criteria' the maximum number of elements in the transfbmied 
vector is five. The present embodiment of the Invention sets the number of elements In the 
output transformed vector to four. 



10 



15 



The data reduction is perfomried In two steps. Refenlng now to Figure 3; all twenty-five of the 
target files 25 (of which four are shown) are processed by an LDA routine 101 that is executed 
by the computer 20 to calculate the LDA transfonnation function gQ. The LDA routine 101 
extracts the -fine' labels from the file names of the target file 25. The firie label is based on the 
first two letters In the name of target file 25. TTie computed gO transform is stored in the 
computer 20. 



20 



Referring now to Figure 4. target files 25 are presented individually to an LDA conversion 
routine 102 that is mn on computer 20. The target file 25 that is input to routine 102 is a two- 
dimensional matrix of 800 rows of Doppler ftame. where each Doppler frame is a 128 element 
vector. Using the LDA transfomiation function gO calculated by routine 101 each row of the 
800 row matrix is a vector with just four elements. Thus the original 128 element Doppler 
frame vector has been reduced to a feature vector with just four elements or features. The 
output of routine 102 is another two-dimensional matrix which has 800 rows but only 4 
columns per row. The output of routine 102 is stored as modified target file 27. The modified 
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target file 27 now represents a sequence of 800 feature yectors for the co-operative example 
targets where each feature vector has just 4 feature (or elements)/ Modified target flies for the 
other 24 co-operative example targets are generated in a iil<e manner using the 
con-esponding target files 25, A total of 25 modified target files are therefore stored In the 
computer 20, each containing data representing the sequence of feature vector of an example 
so-operative target. 



Refening now to Figure 5, modified target files 27 are grouped according to the target group 
fl^om which they are derived. A hidden Markov model (HMM) parameter estimation routine 104 
receives data sets input in turn fl-om groups of modified target files conBsponding to 
personnel, wheeled vehicles and tracked vehicles. It also receives an appropriate HMM 
topology 106, and produces three conresponding sets of HMM parameters at 108. 

Figure 6 shows, a simple HMM of the kind used in speech recognition: it is possible to use It in 
the present invention but better performance is obtained witii the topology shown in Figure 7. 
Figure 7 schematically lllustigtes ttie Hidden Markov Model (HMM) topology implemented by 
the HMM paramet^r estimation routine 104. The HMM topology has Oiree states, 1, 2, 3. 
Conesponding to each state tiiere is a state output distribution which is a multivariate four- 
dimensional Gaussian probability distribution in which ttie probability of observing a given four- 
element vector is given by 



where x is a 4-component vector of a modified data file, //• is a 4-oomponent mean vector 
and E is a 4-by-4 covari'ance matrix. Each HMM state has four mean parameters and 16 
variance parameters (of which only ten are independent as the covariance matrix is symmetric 
by definition). I^or each target group, each of these parameters is calculated from modified 
target files by the HMM parameter estimation routine 104 which implements a HMM estimation 
algorithm of the kind known in the art of speech recognition. Only a fraction of the 800 frames 
from a single modified target file can be used to estimate the parameters of any given HMM 
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State. Thus the 800 frames from a single target file are effectively divided into sub-groups 
corresponding to each state and the parameters. 

In the present embodiment the models are initialised using a simple linear division of the time 
axis into equal segments to give the groups of frames. Other initialisation strategies include 
allocating frames to states on a cyclic basis, with the period of the cycle chosen to optimise a 
goodness-of-fit criterion and/or labelling parts of a database by hand to indicate which groups 
comespond to which frames. 

Subsequently, the well-known HMM re-estimation algorithms perform an autorhatic alignment 
which gives a weighting for the conrespondence of each frame to each state, and weighted 
sums of the feature vectors are used to recompute state mean and variance parameters for 
each state. 

Apart from the state output distributions, an HMM also incorpo^rates a "state transition matrix", 
the non-zero elements of which both define allowed state transitions (the "topology of the 
model) and specify probabilities for occunrence of each allowed transition. Zero elements 
prohibit transitions. Thus the HMM topology may be represented by a diagram such as Figure 
7, in which there is a filled circle con^esponding to each state and an arrow conresponding to 
each non-zero element in the transition matrix. 

The state transition matrix is required to be "row-stochastic", that is the sum of the numbers in 
each row should be unity. It is typically initialised by setting all the non-zero elements In eadr 
row to the same value, a value chosen to satisfy the row-stochastic constraint An alternative 
initialisation for the present invention would be to set the probabilities to small values for "self-, 
loop" transitions (the transition from a state to itselO and larger values for the other allowed 
transitions, thus promoting the desired behaviour of cycling through the states to explain the 
observed data. 

This transition matrix is also optimised by the HMM re-estimation algorithms. In normal use 
those algorithms do not allow zero entries in the transition matrix to become non-zero. .. 
Therefore the initial topology diagram remains valid after re-estimation, although the 
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probabilities associated witli each of tlie transitions will hiave clianged and, in some 
circumstances, some of ttiose probabilities may liave become zero. 

The topology represented in Figure 6 is typical "left-td-righf topology commonly used in 
speech recognition. This can be extended to allow repeated cycling through the sequence of 
states and entry to and exit from any state as illustrated in Rgure 7, which defines the topology 
for three state models used in the example of the Invention described here. 

The topology of Figure 7 has a drawback in that It does not ensure that the full cycle of states 
is observed In recognising a given target a more complex topology is illustrated schematically 
in Rgure 8 which achieves this using the wejUoiown principle of tying of state * output 
distributions. In Figure 8, the states, labelled '1' share a single common probability distribution 
01 and the re-estimation process is modified to pool the con-esponding data; similariy, those 
labelled '2' share another such distribution 02 and those labelled '3' share a further distribution 
03. 

The task of the HMM parameter estimation routine 140 is to optimise the values of the 
transition matrix probabilities and flie state output distributions for each state and each model 
given the tirainirig data. The routine used is well-known in a different field, that of speech 
recognition: it is described for example in "Speech Synthesis and Recognition" by J. N. 
Holmes, Van Nostrand Relnhold (UK) 1988. For each target group, the output of the routine 
140 is a Hidden Maricov IVIodel comprising a state output distiibution parameters for each of 
the states of the model and a state transition matrix. 

Following training as described above, tiie system 10. may be used to recognise and classify 
unknown targets when such targets move wittiin the main radiation lobe of the antenna 12. 
Target recognition Is canrled out as follows. The system 10 is configured witti switch 33 closed 
and switches 34 and 35 open. For a given unknown target, a radar return is digitised to create 
a radar data file as previously described. The radar data file is passed directly to the computer 
20. The routines 52, 53 and 54 previously described are canted out on the radar data file to 
create a target data file corresponding to the unknown target 
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The target file is then operated by routine 102, where the LDA transformation function g() 
calculated during the training phase is used to transfomi each 128 element vector per row of 
the target file 25 to a vector with just 4 elements per row. This then fomis the modified 4- 
element target file which contains 800 frames for the unl<nown target. The modified target file " 
is then processed by the HMM classification routine 206. The present implementation uses the 
well-known HMM Viterbi recognition algorithm described for example in the book "Speech 
Synthesis and Recognition" by J. N. Holmes, Van Nostrand Reinhold (UK) . For a given for a 
given data file, this algorithm efficiently computes and compares matching scores for different 
hypotheses as to the file contents. In the present implementation may be perfomied in "forced 
choice" or "unforced choice" modes. . 

"Forced choice" mode compares the following set of hypotheses for three target classes P, T 
and W: 

{P, T, W, _P, _T, _W, P^ T^ _P^ _T_. _WJ 

where P, T, and W are the target class models and _ represents a background model (no 
target present). "Unforced choice adds a single further hypothesis which is the _ or "no target" 
model. The best scoring hypothesis is reported by the recognition algorithm; the classification 
of the target file is given by stripping the "no target" symbols from the recognition output to 
leave one of the three target symbols (or possibly nothing in the "unforced" case). 

This completes the classification stage where an unknown target is classified into one of the. 
three target classes, i.e. the target is classified as a person, a wheeled vehicle, tracked vehicle 
or as belonging to an unknown dass. 

Refemng now to Figure 9, there is shown an alternative target recognition system of the 
invention indicated generally by 310. The system 310 comprises a decode module 337 for 
decoding data stored in the data recorder 318 and digital signal processors (DSPs) 331, 333, 
333A, 333B and 335. The DSP 331 extracts target Doppler signatures for a single range bin 
as previously described. DSPs 333, 333A and 333B perform Doppler normalisation, 
conversion to Doppler frames and LDA reduction also as described eariier. DSP 335 carries 
out HMM parameter estimation and target classification as previously described. The system 
310 further comprises a display 328. an A-to-D converter 316 and a radar device 314. The 
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computer 332 controls the display 328, the A-to-D converter 316, the radar device 314 and the 
DSPs 331, 333, 333A. 333B and 335 by means of control signals conveyed over control lines 
336. The systems 310 and 10 operate similarly, but in the system 310 processing of radar 
data is earned out by the DSPs 331. 333, 333A, 333B and 335 and not by the computer 332. 

The Invention may use bistatic radar with circular radiation polarisation or any combination of 
horizontal, vertical and circular polarisations. 

The radar devTce may be continuous-wave (CW). The radar device's output power may be 
other than that described above, but should be large enough to obtain a signai-to-noise ratio 
from the target at a desired range sufficient to give satisfactory classification performance. 
The radar device may have a pulse repetition frequency other than 4.096 kHz: however the 
pulse repetition frequency should be sufficiently large to resolve various Doppler frequendes 
assodated with the target and to avoid ambiguous Doppler returns. 

Target recognition systems of the invention may be ground-based, arr-bome or space-borne. 
In the case of embodiments of the invention which move relative to the earth, clutter reduction 
techniques are preferably employed to reduce dutter retums. A scanning antenna may be 
employed. 

Digitised radar retum signals finom co-operative and unknown targets may be stored in a data 
recorder and training and target dassification may be perfomried at a time subsequent to 
collection of training data and radar retums from unknown targets. The transmission frequency 
of the radar 314 may be any frequency used for surveillance radars, e.g. from a few MHz to 
100 GHz. The Invention may use FFT window sizes and window overiaps other than those 
described. Likewise a transfomiation other than LDA could be used to reduce vector size of 
input frames. 

The invention may be set up to dassify more than three classes of targets: a fourth class of 
target e.g. aircraft could be introduced. This would require the appropriate amount of training 
data from the extra dass and estimating the relevant HMM parameters that would define data 
from that dass. The number of states used for each target may be chosen to be different for 
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dififerent target classes. The output distribution of each state could be modelled as a Gaussian 
mixture distribution instead of a simple multivariate Gaussian. 
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CLAIMS 

1. A method of recognising a radar target comprising receiving radar returns from 
a scene and processing tfie returns to produce a Doppler spectmrn characterised in 

5 that the method also includes processing the Doppler spectrum to produce a Doppler 
feature vector and using hidden Markov modelling (HMM) to identify the Doppler 
feature vector as indicating a member of a particular class of targets. 

2. A method of recognising a radar target according to Claim 1 characterised in 
10 that the method includes processing radar returns to obtain a sequence of Doppler 

spectra for each target and producing therefrom a sequence of Doppler feature 
vectors, and using HMM to Jdentlfy the sequence of Doppler feature vectors as 
indicating a member of a particular class of targets. 

15 3. A method of recognising a radar target according to Claim 2 characterised in 
that the method includes using HMM to identify the sequence of Doppler feature 
vectors by assigning to each feature vector an occurrence probability by selecting a 
probability distribution or state from a set thereof associated with a dass of targets, 
multiplying the occurrence probabilities together to obtain an overall probability, 

20 repeating for other probability distributions in the set to detemnine a combination of 
probability distributions giving highest overall probability for that class of targets, then 
repeating for at least one other class of targets and selecting the target class as being 
that which yields the highest overall occun-ence probability. 

25 4. A method of recognising a radar target according to Claim 3 characterised in 
that it includes selecting probability distributions for successive feature vectors on the 
basis of some transitions between distributions allocated to successive feature vectors 
being allowed for the class of targets and others being forbidden. 

30 5. A method of recognising a radar target according to Claim 4 characterised in 
that it includes a preliminary HMM training procedure in which parameters for the 
states or probability distributions and transition probabilities are produced by deriving 
Doppler feature vectors for training data obtained from known classes of targets and 
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calculating thp mean and variance of vectors corresponding to like targets or target 
classes. 

6. A method of recognising a radar target according to Claim 5 characterised in 
5 that the HMM training procedure includes a plurality of cycles through state 

sequences. 

7. A method of recognising a radar target according to Claim 1 characterised in 
that targets encompassed within a single radar range cell. 

10 

8. Apparatus for recognising a radar target comprising a radar receiver for 
receiving radar returns from a scene and processing means for processing the returns 
to produce a Doppler spectrum characterised In that the processing means also 
includes means for processing the Doppler spectrum to produce a Doppler feature 

15 vector and hidden Maricov modelling (HMM) means for identifying the Doppler feature 
vector as indicating a member of a particular class of targets. 

9. Apparatus for recognising a radar target according to Claim 8 characterised in 
that the processing means is arranged to process radar returns to obtain a sequence 

20 of Doppler spectra for each target and to produce therefrom a sequence of Doppler 
feature vectors, and the HMM means is arranged to identify the sequence of Doppler 
feature vectors as indicating a member of a particular class of targets. 

10. Apparatus for recognising a radar target according to Claim 9 characterised in 
25 that the HMM means is arranged to identify the sequence of Doppler feature vectors 

by assigning to each feature vector an occurrence probability by selecting a probability 
distribution or state from a set thereof associated with a class of targets, multiplying the 
occurrence probjabilities together to obtain an overall probability, repeating for other 
probability distributions in the set to detennine a combination of probability distributions 
30 giving higiiest overall probability for that dass of targets, then repeating for at least one 
other class of targets and selecting the target class as being that which yields the 
highest overall occun^ence probability. 
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1 1 . Apparatus for recognising a radar target according to Claim 1 0 characterised in 
that the HMM means is an^nged to select probability distributions for successive 
feature vectors on the basis of some transitions between distributions allocated to 
successive feature vectors ^being allowed for the class of targets and others being 

5 forbidden. 

12. Apparatus for recognising a radar target according to Claim 1 1 characterised in/ 
that the HMM means is arranged to undergo a preliminary training procedure in which 
parameters for the states or probability distributions and transition probabilities are 

10 produced by deriving Doppler feature vectors for fraining data obtained from known 
classes of targets and calculating the mean and variance of vectors corresponding to 
like targets or target classes. - 

13. Apparatus for recognising a radar target according to Claim 1 2 characterised in 
15 that the HMM training procedure includes a plurality of cycles through state 

sequences. 

14. Apparatus for recognising a radar target according to Claim 8 characterised in 
that it is an-anged to encompass targets within a single range • cell. 



wo 02/08783 




PCT/GBOl/03279 



2/7 



Fig.2. 



Decode 
SI 




r 


Extract 
52 




r 


Normalize 
53 


' 


r 


Con 
5 


ivert 
4 



Target file 
25 



SUBSTITUTE SHEET (RULE 26) 




wo 02/08783 PCT/GBOl/03279 

3/7 




SUBSTITUTE SHEET (RULE 26) 



wo 02/08783 



PCT/GBOl/03279 



4/7 



Fig.4. 




Target file 



Modified 
target file 



25 





1 


r 




LDA 
Conversion 


) — ► 








^ 1 





•102 



-^7 



SUBSTITUTE SHEET (RULE 26) 



PCT/GBOl/03279 



5/7 



Fig.5. 



27 










27 




1 







HMM 
Parameter 
estimation 




106 




SUBSTITUTE SHEET (RULE 26) 



wo 02/08783 



PCT/GBOl/03279 



6/7 



Fig.6. 




SUBSTITUTE SHEET (RULE 26) 




SUBSTITUTE SHEET (RULE 26) 



m 




INTERNATIONAL SEARCH REPORT 



ime )nai Application No 

PCT/6B 01/03279 



A. CLASSIFlCA-nONOFSUBJECTMATTBl 

IPC 7 601S7>M1 G01S13/524 



According to Inlamational Palen! Ctas^tlealton (IPC) or to both national classtflcgtlon and IPC 



a RELDS SEARCHED 



Minimum docunnontation saarchod <dafisficaUon systsm toOouwd by dassfficatlon symbols) 

IPC 7 GOIS 



Docurnentatlon saarched other than minlmunn dooimentalion to the extent that such documents em tndidsd In the fields seaiched 



Qectiwitc ddte base consutted during the fadeinational search (name of database and, when practtcat, search terms used) 

EPO-Internal , WPI Data, PAJ, INSPEC 



C. DOCUMB^TS CONSIDEnEO TO BE RELEVANT 



Categoiy Cliallonol document, with IndcaOon, where eppioprtate, of ttie rebvart passages 



HelevanttodelmNa 



NANDAGOPAL ET AL: "Performance of radar 
target recognition schemes using neural 
networks-a comparative study" 
ACOUSTICS, SPEECH, AND SIGNAL PROCESSING, 
1994. ICASSP-94., 1994 IEEE INTERNATIONAL 
CONFERENCE ON ADELAIDE, SA, AUSTRALIA 
19-22 APRIL 1994, NEW YORK, NY, USA, IEEE, 
19 April 1994 (1994-04-19), pages 
II-641-II-644, XP010133782 
ISBN: 0-7803-1775-0 
the whole document 

DE 195 40 722 A (DAIMLER BENZ AEROSPACE 
AG) 7 May 1997 (1997-05-07) 
abstract 

column 1, line 31 -column 2, line 53 
■~ ./- 



1.8 



2,7.9.14 
1,8 

2,7,9,14 



a 



Rutherdooimsnts are listed In the oonlinuafion of box a 



Peienl tamOy mambers are listed In annex. 



* Specbl categories of cited dooinnents : 

'A' docimient dedniig the general state oT the an whicti Is not 

considered to be of parlloular relevanoe 
'E* earfier document but publshed on oranarthe intemaOonal 

nnngdaie 

'L' document which may throw doubts on prtortty daIm(8)or 
wtiich Is cBad to establish the publication date of another 
c&aHon or other special reason (as spedfled) 

'O* document referring to an oral dlsdosun. use. exhibition or 
other means 

'P- docamem published prior to the Intematlanal liDng dale but 
later tlian the prtaifly dale dalmed 



*T* later document published after the (nlemalional filing data 
or prtortiy date and not In conflict with the appllcaibn but 
dtedto understand the principle orthoory underl^hg the 
Invention 

"X* document of particular relevance; the dalmed invention 
cannot be considered novel or cannot be considered to 
Involve an inventive step wtien the document b taken abne 

"Y* document of particular relevance; the dalmed invention 
cannot be considered to Involve an lnvenl^e step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person sMiled 
In the an 

*A* dMumant member or the same patent fandty 



Date of Ihe actual oompletton of the international seardh 



6 November 2001 



Dais of maiing or tne iniemaoonal search lepoit 
16/11/2001 



Name and mailing address of the ISA 

European Patent Office. P3. 5818 Patentiaan 2 
NL-2260HVRil8Vi^iC 
TeL (+31-70) 340-204^ Tx. 31 G51 epo id, 
Fax: (+31-70) 340-3016 



Aulhorizsd ofiioer 



Roost, J 



Form FCT/iaA/zlo (sBoond sheet] [July issz) 



page 1 of 2 



# # 



INTERNATIONAL SEARCH REPORT 


Inie . malAppUcadonNQ 




PCT/GB 01/03279 


C^Contlnuallen) DOCUMEIHS CONSIDERED TO BE RELEVANT 


Category* 


CSation of document. WBh indication,v«tiere appropriate, of the relevani passages 


Rdevani to claim No. 


A 

n 


fiAIiWARD FT AL* "Automatic classification 


1.8 




of environmental noise events by hidden 






Markov models" 






Arnii«;TTr^ ^ppfch and signal processing 






iQQft PRnrrrniNfis of thf iqqa ifff 






TNTFRNATTnUAl CONFFRFNCF ON SFATTLF UA 






IKA MAY NFU YORK NY 






IICA TFFF iJ^ 












ornQ-'^fil ? ypnin?7057Q 






ISBN: 0-7803-4428-6 






abstract * 






dec Lions £ anu o 




A 

n 


^iWLvK n LI nL . IvaOai (.aii^Bt 


1 8 




c> 1 doa 1 T 1 Wf dw 1 UII vi vuuuuei w* 1 a 1 oiiuiaiu 






Trrr TPAN^ArTTOM^ on AFftn^PAOF AND 






Fl FrTRONIC SYSTFMS APRIL 1996 IEEE USA 


















ISSN: 0018-9251 






the Whole document 




A 
n 


MArifFN7TF J 0 FT AL* "Wide-bandwidth 


1.8 




mnhllfi radar for ISAR/SAR radar Imaolna" 






TFF rni 1 noiiTiiM on radar and nicrouavf 






TMARTMfi fOTGFST NO 1994/219) LONDON UK 






MOV 1QQ^ 












1QQd 1 nndnn [K IFF UK 






rIfoH In fho arw^l 1 f^atlon 
v.f 1 tc>u III wMc «Ph 1 luabiuii 






thp whnlp rinnimpnt 




A 

n 


UIIMT FT Al • '*An i nvpctlnatl nn of PI P and 
nuivl LI nL» nil 1 lived t lyab 1 UII ut iLi aiiu 


1 8 




XPIuLUH dwUUdLIu i cpi CdCll UaU 1 Vila OIIU UT 






their potential for combination" 






SPEECH PROCESSING 2, VLSI, UNDERWATER 






SIGNAL PROCESSING. TORONTO, MAY 14 - 17, 






1991, INTERNATIONAL CONFERENCE ON 






ACOUSTICS, SPEECH & SIGNAL PROCESSING. 






ICASSP, NEW YORK, IEEE, US, 






vol. 2 CONF. 16, 






14 Aprn 1991 (1991-04-14), pages 881-884, 






XP010043113 






ISBN: 0-7803-0003-3 






cited In the application 






tne wnoie aocutnenL 





Fbin PCT/lSA/210 (oonniuailon el second daat} (Ju|y isa2) 



page 2 of 2 




INTERNATIONAL SEARCH REPORT 
MwmMlen en patent famUy members 



Inti end Applcatlon No 

PCT/GB 01/03279 



Patant document 


Ribllcation 


Patent t&mlly 


Publication 


cited In search report 


date 


memberts) 


date 



DE 19540722 



07-05-1997 DE 



19540722 Al 



07-05-1997 





Fam PCT/ISMZIO (polamiBiift malt (MV inz) 



